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We study analytically a continuum model for phase-separation in binary polymer blends based on the Flory-Huggins-De 
Gennes free energy, by means of the self-consistent large-n limit approach. The model is solved for values of the parameters 
corresponding to the weak and strong segregation limits. For deep quenches we identify a complex structure of intermediate 
regimes and crossovers characterized by the existence of a time domain such that phase separation is pinned, followed by a 
• preasymptotic regime which in the scalar case corresponds to surface diffusion. The duration of the pinning is analytically 
computed and diverges in the strong segregation limit. Eventually a late stage dynamics sets in, described by scaling laws and 
exponents analogous to those of the corresponding small molecule systems. PACS: 64.75.+g, 64. 60. My, 61.25.Hq 
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(1) ■ I. INTRODUCTION 

s : 

The kinetics of phase separation has been the subject of considerable effort in recent years Among the most 
investigated systems are binary polymer blends beyond their importance for technological applications, such 

■ systems are extremely interesting also from a more fundamental point of view. Due to the intrinsic connectivity of 
, t \ macromolecules, phase separation in polymeric mixtures takes place on space and time scales more easily accessible 
C$ • experimentally than for small molecule systems. Furthermore, long range interactions along the chains greatly reduce 

the size of the critical region, allowing to disregard critical fluctuations in many situations. On the other hand, 
connectivity gives rise to additional complexity in the system, that can in principle lead to a different behavior with 
respect to the small molecule case. Previous investigations H have shown that the global theoretical picture is 
the same for the asymptotic dynamics of a system in the weak segregation limit (WSL), but still many challenging 
£j ■ problems remain unsolved: for example the origin of the pinning phenomenon observed in experiments for off-critical 
quenches and the theoretical study of the strong segregation limit (SSL), which is hard to attack with the traditional 
numerical methods. 

In this paper we consider the most common theoretical model describing the kinetics of phase separation for a binary 
polymer blend: the Cahn-Hilliard equation with the Flory-Huggins-De Gennes (FHDG) free energy functional. We 
\Q ' propose a generalization of the above mentioned model to the case of an 0(n) vector order parameter field ^> a (x, t), 
with a = 1, n, and study the resulting equations in the largc-n limit. The extrapolation of the model to large-n is a 
widely used technique of statistical mechanics, well suited for Ginzburg-Landau-type models, which has revealed quite 
powerful in the study of several systems Q . The large-n limit allows to deal with the nonlinearities of the model by 
means of a self-consistency prescription which effectively linearizes the equations. The solution of the large-n model 
is often a good approximation to the evolution of physical systems with finite n. 

The application of this method to the equation of motion for a symmetric polymer blend allows one to write down 
closed form equations for the main observables and to study analytically the time evolution of the model for all values 
of the parameters. In this way we can present results not only for the often investigated WSL but also for the much 
'■^J ' less known SSL, where a complex preasymptotic scenario can be identified. 

In particular we find that the very late evolution belongs to the same universality class of the small molecule case 
and the typical domain size grows as t x l z with the dynamical exponent z — 4, as usual in vectorial systems. For 
intermediate times, instead, different paths can be taken. More precisely, while for shallow quenches (WSL) the 
evolution closely reproduces what is known for small molecules, for deep quenches (SSL) additional terms, introduced 
in the equation of motion by the polymeric nature of the system, can play a relevant role. In this case, after the 
linear Cahn-Hilliard instability has occurred, the dynamics becomes exceedingly slow and the non equilibrium blend 
remains in a pinned state over an appreciable time interval. This effect is due to the existence of a nonequilibrium 
partially-ordered state which becomes stable for infinitely deep quenches or infinitely long chains. For quenches of 
finite depth the system remains close to this state for a characteristic time t p which is shown to diverge in the SSL as 
a power law or exponentially as the limits of the parameters are taken in different order. Later on it enters a second 
preasymptotic regime characterized by a slow coarsening of domains whose typical size grows as i 1 / 6 . A crossover 
leads finally to the late stage where domains grow as t 1 ^ 4 . 

These results are derived using the large-n model for phase-ordering in polymer blends, which is an approximation 
of the real scalar model. It has been recognized in the small molecule case that this approximation leads to some 
results different from the properties of the corresponding scalar systems. One of them is trivial: the exponents 
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found using such method are valid for vectorial systems and therefore differ from those valid for scalar systems: the 
exponent z = 4 for the growth of domains in the large-n model is known to correspond to z = 3 when the order 
parameter is scalar. Another known discrepancy is that within the large-n approximation the structure factor obeys 
the multiscaling symmetry while systems for finite n exhibit scaling symmetry. This is related to the different nature 
of the phase-transition when n is infinite and in particular to the properties of the equilibrium state the large-n 
model evolves towards; although displaying the correct structure factor such state is not truly ordered: the local 
order parameter distribution remains gaussian instead of becoming bimodal as effect of the formation of domains |g| . 
Therefore a multiscaling structure factor in the large-n model must be interpreted as reflecting scaling symmetry in 
real systems. These differences are well known from the small molecule case and can be easily taken into account, 
but the absence of topological defects in the large-n model may make meaningless the extrapolation of the results for 
the large-n model to the description of actual evolution in scalar systems. More precisely one may wonder whether a 
model with no interfaces (the large-n model) can describe with sufficient accuracy the dynamics in presence of a field- 
dependent mobility distinguishing bulk (where ~ 1) from interfaces (where \ip\ — 0). This problem is discussed at 
length throughout the paper. The key point for the answer is a close comparison of our equations with those of the 
large-n model for ordinary Ginzburg-Landau systems. In such systems the connsection between large-n and scalar 
results is well established. This allows the interpretation of all the complex pattern of regimes and crossovers found 
within the large-n model for polymeric systems in terms of concepts like bulk and surface diffusion even if they cannot 
be defined when n = oo. Hence it turns out that the analytical approximate solution preserves most of the essential 
physics of phase-ordering also for systems more complex than small molecule mixtures. 

The paper is organized as follows. In Section || the usual equation describing phase separation for polymers is 
introduced and the large-n limit model is deduced from it. Section III is devoted to the analytical solution of this 
model and to the identification of the different regimes. The findings of Section |l| are verified and completed by 
means of a numerical solution of the large-n limit equations in Section IV. Finally in Section [v| the results are 
discussed and conclusions are drawn. 



II. THE MODEL 

Model blends are typically described by the Flory-Huggins-De Gennes free energy functional 
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where ^(x, t) is the order parameter field, k B is the Boltzmann constant, a a and a B are the Kuhn lengths of the two 
species, A is an effective interaction distance between monomers, and the Flory-Huggins free energy is [B 
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where Na and N B are the degrees of polymerization of chains A and B; x measures the strength of repulsion between 
unlike molecules and is inversely proportional to temperature. In the following we will always consider for simplicity a 
symmetric blend, for which Na = N B = N and a a = <r B = a. In this case the critical value x^V = 2 separates stable 
states of the blend (xN < 2) from the thermodynamically unstable region (x./V > 2) where the mixture decays in two 
separate phases. The order parameter is related to the volume fraction of A-molecules by ijj(x., t) = 20(x, t) — 1. 

The theoretical description of spinodal decomposition in binary blends is based on the Cahn-Hilliard equation for 
the time evolution of the order parameter field, originally introduced for small molecule systems [Bj 
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where M(ip) is the mobility and n(x, t) is a gaussian white noise with zero average and variance proportional to the 
temperature. For a symmetric blend the mobility 
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has been proposed || where D is a self-diffusion coefficient. With these positions the Langevin evolution equation of 
the order parameter field is given by 
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where we have neglected the thermal noise term since it is possible to show, at least in the analog of this equation 
for small molecule systems ]lO| , pj| , that the temperature is asymptotically an irrelevant parameter below the order- 
disorder line. The scalar model introduced so far can be generalized to the case of an ro-component vectorial order 
parameter ip(x,t) = {-0 Q (x, i)} with a = 1, . . . ,n. As described in Appendix 1, when the number n of components 
diverges the equation of motion for the single component ip a (x, t) (we drop the index a in the following) reads 
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where by definition S(t) =< ip 2 (x,t) > and =< [V^(x, t)] 2 >. 

In the following the dynamics of systems undergoing critical quenches (< ip(x-, t) >— 0) will be studied. In this case 
the quantities S(t) and S*2(i) can be computed self-consistently using the structure factor C(k, t) =< ^(k, t)ip(— k, t) >, 
the Fourier transform of the real space pair connected correlation function, through 
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where d is the spatial dimension of the system and q is a phenomenological ultraviolet momentum cutoff. The evolution 
equation of the structure factor can be obtained from Eq. (0) as 
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Eq. (^0() together with the self-consistency relations (|8|) and (|9j) constitute the integro-differential equations governing 
the dynamics of the FHDG model in the large-n limit, and will be the object of our study. Some differences occur 
between Eq. ( |l0| ) and its analog for small molecule systems [||: the overall time-dependent factor 1 — S(t) reflects the 
presence of a field-dependent mobility as opposed to the constant value usually taken; the terms proportional to a 2 
are a consequence of the order parameter dependence of and do not appear in the equation for small molecules; 
the logarithmic form of /fh is different from the usual Ginzburg-Landau quartic potential. 

As we will see, in phase-separation of binary polymer blends a key role is played by the product xJV, setting how 
strongly segregated the two species are. When xN ^ 2 the absolute value of the equilibrium order parameter in the 
coexisting phases is much smaller than 1, indicating that in A-rich regions a high concentration of B molecules is 
present and vice versa. This is the Weak Segregation Limit (WSL), which is also called shallow quench condition, 
since it is usually realized by lowering T just below T c . For deep quenches instead xN ^S> 2 (Strong Segregation Limit, 
SSL); the order parameter equilibrium value is very close to ±1 because separated domains are almost pure. 
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III. ANALYTICAL STUDY OF THE MODEL 



In the following a high temperature disordered initial condition C(k, 0) = A will be considered. Upon introducing 
the three quantities 

L(t) = (\NDa 2 t\ 4 (11) 
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Eq. ( |l(]| ) can be formally integrated, yielding 

C(k,t) = Aexp {k 2 [A 4 (t) + L\t)] (2k 2 m (t) - k 2 )} 
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is the position of the peak provided L 2 > 0, as is the case for sufficiently long times in the phase-ordering region. 

The introduction of the three quantities (0), (|l2]) and ( [l3"|) allows the description of the dynamical evolution in 
terms of their competition. Depending on their relative size, the system exhibits different properties in subsequent 
time regimes, as will be described in detail below. 



A. Early stage 



For short times after the quench, assuming the initial fluctuations A are not large, S(t) and S*2(i) can be neglected 
in Eq. (10) and the system exhibits the usual linear behavior of phase-ordering. The time evolution of A(t) and £(t) 
is easily computed 
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As a consequence, in complete analogy with linear behavior in small molecule systems, the position of the maximum 
remains constant 



while its height grows exponentially fast 
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During the very early stages of this linear regime (for times t -C t* , where t* will be determined below) the structure 



factor obeys an approximate scaling form. This can be shown by considering that, according to Eq. (14), C(k, t) 
decays for |k| > k m over the typical distance W(t) = [A 4 (t) + L 4 ^)}- 1 / 4 = [ND( X \ 2 + a 2 /%)t]- x / 4 . Since for small 
times W(t) is very large, the integrals defining S(t) and Si(t) are dominated by the contributions for large momenta. 
For large k, k 2 n can be neglected in Eq. ( |l4| ) and the structure factor can be written as 

C(M) = Aexp{-/c 4 [A 4 (t) +L 4 (t)}} = Aexp {-g[kL(t)} 4 } (20) 

with g constant. Hence C(k, t) exhibits for k 3> k m a scaling form with respect to the growing length L(t); S(t) and 
Sz(t) can be computed easily, yielding 

S(t)~AL(t)- d S 2 (t) ~ AL(t)- (d+ V (21) 

consistently with the assumption that they are small. This very early scaling regime is completely analogous to 
the regime found in small molecule systems for very short times before the usual Cahn-Hilliard linear regime [ jT2| . 
Its physical origin is the presence of totally uncorrelated fluctuations in the initial state, creating large gradients 
in concentration between neighboring regions: the square gradient is the dominating contribution to the excess 
free energy. Then the system lowers its free energy by reducing everywhere the local order parameter so that the 
contribution of the square gradient is reduced. In this way the energy associated with the local potential grows, but 
as long as it is much smaller than the other contribution it does not affect the evolution: the behavior is diffusive as if 
the system would be at x_/V = 2. This type of evolution ends when the global reduction of the local order parameter 
ends up increasing the total free energy, because the local contribution grows more than the decrease in the square 
gradient one. In terms of the structure factor this happens when the peak position is no longer much smaller than its 
width W(t), i.e. for t = t* such that 

fc ro (r)~ [A*(f) + £ 4 (t*)]~ 1/4 ( 22 ) 

from which 

_ 64jV( X A 2 +a 2 /9) 

- D( xN -2r (23) 

From Eq. ( p3| ) we conclude that the duration of the very early scaling regime diverges for infinitely shallow quenches. 
In the same limit k m vanishes and therefore the scaling form (^oj ) holds down to k = 0. After the crossover time t* 
the usual linear behavior sets in and the order parameter saturates exponentially fast to the local equilibrium, leading 
to an exponential growth of S(t) and S 2 (t). 

B. Pinned regime 

The linear regime continues until a time to such that S(to) is close to the equilibrium value iS(oo) which corresponds 
to the minima of the Flory-Huggins local potential. The subsequent behavior of the system strongly depends on the 
value of x, N and a. In the WSL 5(oo) <C 1 and therefore A(t) = (NDx^H) 1 ' 4 ~ Lit); the system enters immediately 



the asymptotic stage described in Sec. [II D. In the SSL instead, the dynamical evolution of the large-n equations 
enters a quasi-stationary regime which extends over a time domain diverging when X N — > oo and X /u 2 ~~ * °°; this 
case will be referred to as the pinning limit. During this time interval no appreciable evolution of C(k, t) is observed 
and the blend is practically pinned in a configuration characterized by phase-separated domains of finite size. This 
phenomenon reflects the existence of a static non-equilibrium configuration becoming stable in the pinning limit. In 
order to see this let us consider the static solutions of the model. For general values of the parameters it is clear that 

C(k,oo) = (2 7 r) d 5 oo( 5(k) (24) 

is a static solution (dCJdt — 0) of Eq. (|l0|), where the choice of = S(oo) as prefactor is dictated by the self- 
consistency condition (|8|). The value of 6*00 is fixed by requiring the solution (|2J) to be a minimum of the FHDG free 
energy. The resulting condition is 

1 -lnf 1± #Vx = (25) 



NSl V 1 - 5| 
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which in the scalar case indicates that the order parameter lies on the minima of the local potential. This is the 
equilibrium state towards which the system evolves for all finite values of xiV and x/ a2 i as can be checked using the 



asymptotic results of Sec. HID. A different situation occurs instead when \N and x/tx both diverge. In this limit 
the equation of motion reads 

xND[l-S{t)]k?l-l + \*Ac{Kt) (26) 



dt WJ 1 4 

Integrating one obtains for all times 



C(fc^) = Aexp|-fc 2 A 4 (t)^ 2 -^| (27) 
Therefore a static solution requires A(i = oo) = Aoo = const, yielding 

C(k,t) = C p (k) = Ae W ^-k 2 At ( fc2 -^)} (28) 

We notice from Eq. (Jl^) that a constant A(t) implies S(t) = 1 and this sets the value of A m via the self-consistency 
condition 

(2^W = 1 ( 29 ) 

Hence in the pinning limit the system does not evolve towards the free energy ground state corresponding to 
complete order, as revealed by the lack of the Bragg peak at k = in C p (k). Instead a state partially ordered over 
a typical length A:" 1 = 2^/2X is dynamically generated; despite having a higher free energy than the equilibrium 



configuration, the state (28129) is strictly asymptotic in the pinning limit. When \N and xl® 2 are large but finite, 



the additional terms in Eq. ([TOP destabilize the pinned state: the system gets trapped around it for a time domain 
diverging in the pinning limit. The duration r p of the interval during which no coarsening practically occurs can be 
easily estimated, as reported in Appendix 2. It turns out that this time strongly depends on how the pinning limit is 
approached, i.e. on the order of the limits \N — > oo and x/ a2 ~^ 00 ■ If one takes x/c 2 — > oo first 

T p ~exp{a X N} (30) 

where a is a known constant, while when — > oo first 

1/2 



r P ~ (£) (31) 



This twofold behavior reflects the different terms that can destabilize the pinned state. By comparing Eq. ( |To| ) and 
Eq. ( p6] ) it turns out that Eq. ( JToj ) contains two kinds of additional terms: the logarithmic contribution and the 
terms proportional to er 2 . The former is proportional to 1 /N and hence destroys pinning when xN is finite leading to 
Eq. ( p0|); the latter are active when a is nonvanishing and cause the decay of the pinned state over the characteristic 
time (pip . 

It is important to stress that the pinning regime occurs when phase separation has already taken place but is still 
incomplete and, therefore, has nothing to do with the usual metastability present during phase ordering, which decays 
via nucleation. 



C. Surface diffusion regime 

In the SSL, after the end of pinning the structure factor has the form (|l4|). During this regime C(k,t) is sharply 
peaked around k m (t); this allows the evaluation of S(t) and Sa(£) by saddle point technique, yielding 

S(t) - /T d (f)e^m+A4( t )] ( 32 ) 

and 

S 2 (t) ~ £- 2 {t)S(t) (33) 
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Requiring the saturation of 5(f) to its equilibrium value 5oo and neglecting L(t) with respect to A(t), for sufficiently 
large x/ over(j2 j one nas 

C d (t) ~ e" 1 ^) (34) 
According to its definition (p"3|), to leading order in — 5(f), C 2 (t) can be written as 

[ 5 1 / 2 (r)-5^ c /2 ][l-5(r)] a 2 5 2 (r) \ 

' 45V^ "9" [1 - 5(r)] 2 J dT (d5) 

with 6 = 2/[N(l — Soo)]. Then, neglecting for sufficiently small Kuhn length a the term containing Si(€) in Eq. (|3^), 
and letting 5oo — 1, we obtain 

£ 2 (i) ~ £ 2 (i ) + / [1- 5? (r)] 2 dr (36) 
where t is the crossover time from pinning to this regime. In the same way, one finds 

A 4 (f)~ A 4 (f )+ 55 (r)]dr (37) 

./to 

Assuming that C(t) and A(t) diverge with time, we neglect the constant terms in Eqs. ( j36"l ) and (|37|). Eqs. (|34|) , (|36|) 
and ( |37| ) admit then the solution 

A(t)~t* £(t)~t*(logt)* (38) 

with f — 5 1 / 2 (i) ~ < -1 / 3 , consistently with the assumption of diverging C(t) and A(i). As a consequence, using 

Eq.©,M*)~*- 1/8 0°gt) 1/4 - 

During this regime the dynamics is governed by A(i) which dominates over Lit). With the help of definitions jll| ) 
and (|l^) it is clear that during this stage the system behaves as if the square gradient coefficient k was independent of 
-0, i.e. as if a = 0. With that condition Eq. ( ^0| ) becomes perfectly analogous to the equation of motion for the large-n 
approximation of a Ginzburg-Landau system with a field dependent mobility M(ip) = 1 — ip 2 ; it is easily recognizable 
that also in this case the analytical solution yields k m (t) ~ f -1 / 6 . This similarity with the Ginzburg-Landau system 
helps understanding the physical meaning of this regime. The scalar case for the small molecule Ginzburg-Landau 
problem with non constant mobility has been studied by means of a Lifshitz-Slyozov approach jl3| and numerical 
simulations [^4|. The outcome of such investigations was that when n — 1 the typical length grows as i 1 / 4 and the 
dominating growth mechanism is surface diffusion. This leads to the conclusion that also for a polymeric mixture the 
time regime with k m (t) ~ t~ 1//6 in the large-n model reflects a surface diffusion regime in the corresponding scalar 
model with k m (t) ~ i" 1 / 4 . This is why (with an abuse of language) we termed this regime as surface-diffusive: clearly 
no surfaces exist in the large-n model and no diffusion along them takes place. Nevertheless the behavior of the 
large-n model clearly reflects the prevalence of this growth mechanism in the corresponding scalar system. 



D. Asymptotic regime 

The behavior illustrated in the preceding section is not yet asymptotic, since A(t) ~ t 1 / 6 cannot dominate L(t) ~ t 1 / 4 
forever. When this occurs, the system enters the very late stage of its temporal evolution. What follows is valid also 
for the WSL; in such case L(t) and A(t) are always proportional to i 1 / 4 and the asymptotic regime begins right after 
the linear one. 

The analytic treatment of this late regime is again based on the saddle point technique, yielding Eq. (|3^), and the 
request that the order parameter field lies on the ground state manifold, S(t) = 5oo. From Eq. ( |i"l| ) and Eq. ( O ) 
A(t) - L{t) - t 1 / 4 then Eq. @ gives 

C 4 (t) ~L 4 (t)logt (39) 

and therefore 

k m {t) ~ (log t/t) 1 / 4 (40) 
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Here the two lengths L(t) and k m (t) diverge in the same way up to a logarithmic factor yielding multiscaling 
precisely as for phase-ordering in large-n ordinary mixtures H] and showing that the two models fall into the same 
universality class. Again we can use the comparison with phase-separation in ordinary mixtures to extrapolate results 
for scalar polymer blends: on the basis of what is known Jll],[l5| about small molecule systems, a scaling regime is 
expected to be obeyed also for polymers for finite n, characterized by a single diverging length growing as t 1 ^, with 
z = 4 for vector order parameter and z — 3 in the physically relevant case n = 1. The coarsening mechanism prevailing 
during this stage in scalar systems is bulk diffusion: although strongly suppressed in the SSL, it is not vanishing and, 
being associated with a faster domain growth, finally dominates over surface diffusion. 

It is interesting to remark that the model for polymers has the same asymptotic behavior of small molecules, but for 
non trivial reasons. In ordinary mixtures L(t), the dominating length during the late stage, is formed by the product 
of a constant mobility times a constant square gradient coefficient. For polymer blends L(t) is the result of a non 
constant M(i/j) times the field dependent part of k(i/j): the expression for L(t) is the same of the Ginzburg-Landau case 
because the order parameter dependence in the two factors cancels out. Therefore, even if the asymptotic behavior 
of this model for polymer blends is the same as for small molecules, it is not correct saying that the field dependence 
of M and k is irrelevant. 



IV. NUMERICAL RESULTS 

In this section we present the results of the numerical solution of the large-n model. The solution is performed by 
simple iteration of the discretized version of Eqs. with ^ = 3 and 1024 values of k. The diffusion coefficient D 

is chosen equal to 4 and the number of monomers is fixed to N = 0.25 x 10 5 . The value of the parameter \ is changed 
over many orders of magnitude, so that we can clearly distinguish the different time regimes. 

We start by considering the very early stages. The behavior of S(t) and 5*2 (t) for very early times is shown in 
Fig. [I]. The initial decay of S(t) follows very accurately the power-law t -3 / 4 found analytically. The same agreement 
is found for S2(t) with the decay t~ 5 / 4 . Notice that the crossover time is very close to the estimate based on Eq. (|23|): 
r~4xl0~ 2 . 

With regards to the following stages, in Fig. || 52 (t) and k m (t) are plotted versus time for values of the parameters 
in the weak segregation limit. The linear behavior is clearly visible, characterized by a constant position of the peak. 
It is followed by a sharp transition to the asymptotic regime, during which the two plotted quantities decay as power 
laws. If k m (t) is fitted with (t/ logi) 1 / 2 , the computed exponent is 1/z = 0.253 ± 0.003, in good agreement with the 
theoretical value z = 4. 

For deep quenches in the strong segregation limit, the situation is quite different (Fig. ||). At the end of the linear 
regime, the onset of the very late stage dynamics is preceded by the two preasymptotic behaviors mentioned in the 
previous section. During the first one the system undergoes an almost complete stop; therefore k rn (t) remains at 
the value of the linear regime, while S2{t) shows a plateau which extends over many decades. This temporary stop 
in the evolution of the system is even better illustrated by plotting directly the structure factor C(k,t) for different 
times during the pinning regime (Fig. ^). These curves are compared with the analytic expression ( p8[ ) of C p (k) in 
the pinning limit. Later the dynamics restarts and is dominated by the time dependent mobility. The peak position 
and £2(2) go to zero as power-laws, in good agreement with the expected behavior, k m ~ i" 1 / 6 and S2 ~ t^ 1 ^ 3 . The 
agreement is not perfect because the system is already crossing over to the asymptotic behavior. The onset of this 
last regime can be delayed by making x/ a2 bigger, but this would also increase the duration of the pinned stage, 
making the surface diffusion regime numerically unreachable. Finally, on times longer than those shown in the figure, 
both quantities smoothly cross over to the asymptotic behavior, which is the same of the WSL. 

The duration t p of the pinning is displayed in Fig. M. In the upper part r p is plotted versus x/ <j2 f° r N strictly 
infinite, showing a power law behavior whose measured exponent is 0.49 ± 0.01, in very good agreement with the 
analytical estimate of Eq. In the lower part the same quantity, computed for a = 0, is plotted versus %A, 

displaying an exponential dependence as predicted in Eq. (pp|). 

In the end, all figures confirm the analytical results discussed above and the existence of a complex structure of 
intermediate regimes and crossovers, as summarized in Table 1. 

V. DISCUSSION 

The solution of the large-n model for phase-separating polymer blends leads naturally to a comparison with the 
analogous results for small molecule systems. In this way we can identify which of the modifications introduced by 
the macromolecular nature of the blend components are relevant. We consider the effect of three modifications: 1) 
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the mobility depending on the local order parameter and in particular vanishing in pure phases; 2) the local potential 
having a double well form, but a logarithmic expression, as opposed to the usual polynomial; 3) the square gradient 
coefficient in the free energy having an additional contribution depending on the local order parameter, giving rise to 
two new terms in the chemical potential. The first two differences are actually not restricted to polymer blends and 
can be considered also for small molecules; the third is instead strictly related to the macromolecular nature of the 
mixture. 

In the WSL our results confirm what was already known from previous numerical simulations of the full continuum 
equation The system belongs to the same universality class of small molecule blends and all additional terms 
of the model turn out to be irrelevant during the whole dynamical process. This is not surprising when the scalar 
system is considered: in the WSL separated phases are not pure: the equilibrium order parameter in phase-separated 
domains is far from -1 and 1; in such conditions M(ip) and n(ip) have small local variations negligible compared to 
the average constant terms. 

The situation is much more interesting in the SSL. It turns out that all three modifications are relevant in this 
case. The logarithmic expression for fpn pushes the minima of the local potential close to -1 and +1 exponentially 
with xN • For these limit values the mobility vanishes, and the evolution is pinned. During the subsequent regime 
the non constant mobility is relevant since the evolution is governed by A(t) which owes its time dependence to the 
order parameter dependence of M(tp). Finally, as already pointed out in Sec. |lll D , the asymptotic stage is governed 
by the growing length L(t) which is formally the same of the small molecule case but is actually the result of the field 
dependence of M(^) and K,(ip). 

Our aim is also to make statements about the real systems, not only about the large-n approximation to their 
temporal evolution. Therefore a word must be said about the delicate problem of the connection between the systems 
we want to study (scalar order parameter) and those we are able to solve analytically (vectorial order parameter with 
an infinite number of components) . For the small molecule case we already know that some properties of the solution 
for large-n model do not hold for scalar systems. One of them is the multiscaling symmetry of the structure factor for 
long times: for finite n scaling symmetry holds. Another difference is the value of the dynamical exponent z which is 
known to be 3 for scalar systems while is 4 when the order parameter is vectorial (including the large-n model). 

These differences are known and can therefore be easily taken into account. More dangerous may in principle 
be another difference between the large-n and the corresponding scalar model: the latter forms ordered domains 
separated by well defined interfaces; the former does not support interfaces and actually evolves towards a state that 
is not truly ordered ||. This difference could be critical for a polymer blend in which interfaces play a key role, 
through the field-dependence of the mobility and of the square gradient coefficient. Nevertheless we believe that in 
this case the picture provided by the large-n model is a close representation of what actually goes on in scalar systems. 

This conclusion relies on the comparison of our results with those obtained for small molecule Ginzburg-Landau 
systems within the large-n approximation with constant and non constant mobility. All the behaviors we find can 
be found also in large-n models for Ginzburg-Landau systems, where they are interpreted as the result of different 
physical mechanisms governing growth. By analogy we can describe our large-n results as the effect of the interplay of 
competing coarsening mechanisms for scalar polymer blends. The pinning regime, the subsequent regime characterized 
by z — 6 and the asymptotic stage are all very clearly interpreted in terms of the growth processes occurring in scalar 
polymer mixtures. 

In particular, two are the mechanisms driving coarsening in binary blends. The first is bulk diffusion, also said 
Lifshitz-Slyozov or evaporation-condensation mechanism: A-molecules evaporate from high curvature regions of A-rich 
domain interfaces; they diffuse in B-rich regions and condensate on A-rich domains with lower curvature. This process 
makes smaller domains shrink and larger grow; it is associated with a t 1 ^ 3 growth law becoming i 1 / 4 for n > 1. The 
competing mechanism is the diffusion of molecules along domain surfaces in order to minimize the interfacial energy. 
This process has the effect of changing the shape (but not the volume) of single domains and is associated with z = 4 
(z = 6 in the corresponding large-n limit). The slower growth law explains why surface diffusion is not observed in 
usual small molecule systems and in polymers in the WSL: bulk diffusion always prevails. For deep quenches instead 
both mechanisms are slowed down, but in different fashions. Surface diffusion depends little on temperature and is 
only weakly suppressed when T — * 0. On the other hand the evaporation process needed for the Lifshitz-Slyozov 
mechanism is activated and therefore exponentially inhibited for deep quenches: its probability is proportional to 
exp(— AF/kBT) and the free energy change involved by the evaporation of a macromolecule is AF ~ ksTxN. In this 
way we can interpreted the succession of stages occurring in the SSL: when both growth mechanisms are inhibited 
the system is pinned; later growth starts, driven by surface diffusion, which is slow but only weakly suppressed. 
Eventually bulk diffusion prevails and phase-separation enters its late stage. 

Within this context it is difficult to understand why the pinning limit requires xl — > oo in addition to %iV — * oo. 
It is plausible that the additional condition is required only for the large-n model and does not apply to scalar systems. 
This is suggested by the observation that the condition x/cr 2 — > oo is needed in order to neglect the term proportional 
to £2(4) hi Eq. jlO|). Such term appears in the equation with same role of those derived from the local potential in 
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the free energy. However it actually comes from the nonlocal part of the free energy: it becomes "local" (i.e. not 
proportional to k 2 ) only as effect of the large-n limit. It is very likely that in the scalar case the evolution freezes 
even if x/ " 2 < 00 ■ 

We finally discuss the relevance of the previous results in the interpretation of the pinning phenomenon that has 
been observed in off-critical quenches. Experiments show that some polymeric mixtures quenched in the unstable 
region of their phase diagram dramatically change their behavior depending on the average concentration of the blend 
components |l(|[l7|. When the concentration is critical growth proceeds as usual. When concentration is sufficiently 
off-critical coarsening starts but later stops, before the system reaches equilibrium, in a frozen configuration with 
partially separated phases. The specific mechanism responsible for this phenomenon is still poorly understood and 
this topic has been the subject of discussion recently |jl^,^9|. The conjecture that inhibition of bulk diffusion due 
to free energy barriers may play a fundamental role has been put forward Jl6| ] but so far no convincing test of this 
hypothesis has been done: direct numerical integration of the full equation of motion is easily performed only in the 
WSL and no pinning has been detected |t9|; for deep quenches spurious numerical instabilities arise. Using the large-n 
limit approximation we are able to investigate the strong segregation limit. From our study a plausible explanation 
of the experimental evidence comes out. For extremely deep quenches all growth mechanisms are suppressed and 



the system is pinned in a configuration out of equilibrium. This is the pinning described in Sec. Ill B and does not 
depend on concentration, i.e. it happens also for critical quenches. It is very unlikely that this kind of pinning 
is observed in experiments, since it probably requires unrealistically low temperatures. The pinning phenomenon 
observed experimentally in instead more likely related to intermediate values of %-ZV, such that bulk diffusion is 
inhibited while surface diffusion is not. This would explain both the unarrested growth for critical quenches and the 
freezing for off-critical ones. When the concentration is critical an interconnected pattern is formed for both phases and 
surface diffusion can drive the system to macroscopical phase-separation. When the quench is sufficiently off-critical 
instead, the minority phase forms non percolating droplets embedded in a matrix of the majority phase. Surface 
diffusion can only lead to a partial phase-separation and coarsening stops when droplets are spherical. However, on 
much longer times, the residual bulk diffusion should drive the system to complete phase-separation. In order to 
confirm this scenario further work on the numerical solution of the scalar order parameter equation in the SSL is in 
progress. More experiments, aimed at verifying the prediction that coarsening should restart for very long times after 
the pinning in off-critical quenches, would also be very helpful. 
We thank Marco Zannetti for an interesting discussion. 

APPENDIX 1 

We consider the Flory-Huggins-De Gennes free energy functional F(ip) and the mobility M(ip); in order to generalize 
them to the vector order parameter case we require FV(V') an d Mv(ip), the vectorial counterparts of the free energy 
and of the mobility, to be O(n) symmetric functions of ?/>(x, t). With this position the field dependence occurs through 

the modulus of vector quantities, namely \tp (x, t) | = YHp=i ^( x i f ) an< i |V^(x, t)| = < ^^ =1 [V^Va( x J *)] 2 f 
In the large-n limit if one requires the single component ^ a (x, t) to remain finite, the square modulus of vector 
quantities must be normalized by 1/n, in order to keep it finite. Hence the whole field dependence of i*V(V0 an d 
My(ip) occurs in the vectorial case through |i/>(x, i)| and |V?/>(x, where ijj(x,t) = n~ 1//2 ^(x, t). Moreover one 
requires Fy{ijj) to be an extensive quantity in the number of components n. In summary, a proper generalization to 
the vector case is achieved by substituting everywhere ip(n,t) and |V^(x,t)| with |^>(x, i) and |V^>(x,t)| respectively 
inEqs. (0), (§, @, (|), and multiplying F(V>) by n. We obtain 

Mfi^r^tfMm) , Kimm{lc , t) A (4i) 



k B T J \ k B T 

where k(IV'I) an d /fh (li^l) are still given by expressions @ and @. The Cahn-Hilliard equation for the time evolution 
of a generic component a of the vector field ip(n, t) reads 



dt 



blp a 



(42) 



where My(-0) = M(|-0|) is given by Eq. (||) and we have neglected thermal noise as discussed in Sec. 0. Then, 
considering a symmetric blend, the Langevin evolution of i() a (x.,t) is obtained 
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dip a (x,t) ND 



dt 



■V-Ul-^(x,*)| 2 )V 



j_ / i + |^(x,t)| \ M*,t) 

AN n Vl-|?(x,t)|J |^(x,i)| 



~V>a(x, t) ~ U'x + = 



V 2 ^ a (x,t) + 



Q-Va(x,t) 

9(l-hKx,*)| 2 ) 2 



|VV a (x,i)| 5 



(43) 



For n = 1 one recovers Eq. (|6j). In the large- n limit, summing over vector components averages the system over an 
ensemble of configurations, and hence 



— 1 n 

lim |V(x,t)| 2 = lim - VV>|(x,i) =< Va(x,*) >=S(t) 



(44) 



0=1 



where < ... > denotes the ensemble average, translational invariance has been assumed and S(t) does not depend on 
a due to internal symmetry. Analogously 



_ 1 " 

lim |V^ Q (x,t)| 2 = lim - V[V^a(x,i)] 2 =< [V^(x,t)f >= S 2 (t) 

n — >oo n— »oo fl z — ' 



(45) 



/9=1 



Hence Fourier transforming to reciprocal space, the evolution equation for the order parameter field reads 



™[i _ W ,„ . ^ 



In 



, ° 2 &(*) 
9 [1-5W] 2 

where the component index a has been dropped. 



9[1 - 5(t)] 



fc 2 U(k,t) 



(46) 



APPENDIX 2 



In this Appendix we derive the expressions (30) and ( |3i| ) for the duration of the pinned stage when the pinning 
limit is approached. We define t p with reference to the behavior of the quantity S^i). As can be seen in Fig. ^ S%(t) 
displays a plateau during the pinned stage. More precisely it reaches a maximum for t = t p at the end of the linear 
regime and decreases extremely slowly until the crossover time t = t p + 6 p , when the pinned stage ends and the system 
enters the subsequent time evolution characterized by a more rapid (power law) decrease of Sa(t). A quantitative 
definition of 8 p can be obtained from the relative variation of S 2 by letting 



S2(tp) — S2\t p + 8 P ) 



(47) 



S 2 (t p ) 

where e is an arbitrarily fixed small number. Since is approximately constant during the pinned stage one has 

dS 2 (t) 



s 2 {t P + P ) = S 2 {t p ) + 



dt 



(48) 



and therefore, using Eq. ( f47| ) 



-eS 2 (t p ) 



' dS 2 {t) 




dt 





(49) 



Op is the actual duration of the pinned regime but it should be noticed that in the pinning limit dC/dt is proportional 
to X-W) which goes to infinity. All times are divided by this factor and hence vanish. In order to compare the duration 
of the pinned stage for different values of xN , P must be rescaled by the appropriate intrinsic time factor l/\N; we 
therefore define the duration t p of pinning as 
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r P = X N9 p (50) 
The derivative in Eq. (|49|) can be computed by considering that C(k,t = t p ) ~ C p (k) defined in Eq. (|2|) so that 



dS 2 (t) 



Of 



d d k 2 <9C(k, t) 

7 «C 



(27r) d 9< 



-XJV£>(1 - 5 p )(A 2 5 6p - S 4 p/4) 



where 



<Tk 



(27T) 



-^k n C p {k) 



Snp — Sn(tp) - 

are known quantities. By inserting Eqs. ([l9|) and ( pl| ) in Eq. ( pp| ) one obtains an expression for r p 

^ ejjp 

Tp ~ £>(l-S p )(A 2 S 6p -S 4 p/4) 



(51) 



(52) 



(53) 



where only the value S p of S(t) during the pinned stage remains to be determined. This is calculated by imposing 
that 



dS(t) 



at 



o 



(54) 



When evaluating this condition the outcome depends on how the pinning is approached, i. e. on the order of the 
limits %iV — > oo and x/ (j2 ~ ¥ 00 • When one takes x/^ 2 ^ 00 with large but fixed N, one has 



= - X ND{1 - S p ) 
With simple algebra one obtains 



-1/4- 



•In 



5*2p + A Sip 



1-S p = 4expl- X N 



(S2p — 4A 2 S , 4 P ) 



S 



2p 



and therefore 



r p ~ exp <^ x^V 



(<S2p — 4A 2 5 , 4 P ) 



'2/j 



Letting instead xN — > oo with non vanishing cr 2 



= - X JV£>(1 - 5 P ) 



-1/4- 



a 2 5- 



2p 



S 



2p 



yielding 



and therefore 



i-s p = 



x 9(i - s p y 

' ^2p 

9(5 2 p/4- X 2 S 4p ) W 

1/2 



A 2 + 



a 2 1 



X 9(1 - S p ) 



S 



xr 1/J 



(5) 



rp ~ (£) 



(55) 



(56) 



(57) 



(58) 



(59) 



(60) 
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FIG. 2. Top: Log-log plot of S2 (i) vs t for a quench in the weak segregation limit. The values of the parameters are 
N = 0.25 x 10 s , xN = 2.1, A = 1/2, cr = l. Bottom: The same plot for k m (t). 
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FIG. 3. Top: Log-log plot of vs t for a quench in the strong segregation limit. The values of the parameters are 

N = 0.25 x 10 5 , xN = 10 12 , A = 1/2, a = 1. Bottom: The same plot for k m (t). 
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FIG. 4. Plot of C(k, i) vs t for the same parameters of Fig. (j3|) and two different times separated by two decades, compared 
with the analytical expression (Eq. psh). 
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FIG. 5. Top: Log-log plot of t p vs x/°" 2 f° r N = oo showing the power law divergence of the pinning duration. Bottom: 
Linear-log plot of r p vs \N for a — 0, displaying that in this limit t v diverges exponentially. In all cases t p was determined by 
choosing e = 10 -4 . 
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TABLE 



Regime 


Early scaling 


Linear 


Pinned 


Surface Diff. 


Asymptotic 


S(t) 


t -d/i 


exp(i) 


const 


const 


const 


s*(t) 


t -(d+2)/4 


exp(t) 


const 


t- 1 /3(logt)-V2 


t- 1 / 2 (logi)- 1 /2 




const 


const 


const 




(i/logi) 1/4 


C(t) 


t l/2 


t l/2 


const 


tV6(lo gt )l/4 


(tlogt) 1 ^ 


A(t) 


t l/4 


t l/4 


const 


t l/6 


t l/4 



Table 1. Summary of the time dependence of the important quantities during the different stages. The third and 
fourth time regimes are observable only when \N 3> 1 and xl ' j2 3> 1. 



1G 



